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Abstract 

In many vibration minimization control problems for high rise buildings subject to strong 
earthquake loads, the emphasis has been on a combination of minimizing the displacement, 
the velocity and the acceleration of the motion of the building. In most cases, the 
accelerations that are involved are not necessarily large but the change in them (jerk) are 
abrupt. These changes in magnitude or direction are responsible for most building damage 
and also create discomfort like motion sickness for inhabitants of these structures because 
of the element of surprise. We propose a method of minimizing also the jerk which is the 
sudden change in acceleration or the derivative of the acceleration using classical linear 
quadratic optimal controls. This was done through the introduction of a quadratic 
performance index involving the cost due to the jerk; a special change of variable; and 
using the jerk as a control variable. The values of the optimal control are obtained using 
the Riccati equation. 

1. Introduction 

Classical optimal controls in buildings typically consider the control of displacements and 
velocities [Loh and Ma, 1994, Soong and Yang, 1988]. There is minimum consideration 
for the control of acceleration [Yang and Li, 1991] or the rate of change of acceleration, 
otherwise known as jerk IFinney and Thomas, 1994],. This imposes limitation on the 
effective control of motion sickness resulting from jerks. There is therefore a need for 
more precision about the modeling of control forces in buildings since seismic excitation is 
jerky. For a building structure subjected to earthquake excitations, a jerk minimization 
method for vibration control is proposed. 

2. Equations of motion 

We consider a MDOF system of the form 

Mx + Ci + Kx = Du + Ef (1) 

where x is the n-dimensional displacement vector, M, C, and K are respectively the n x n 
mass, damping and stiffness matrices, u is the m-dimensional control vector forces and f is 
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the r-dimensional external excitation vector forces, whereas, D and E are respectively the 
n x m and n x r location matrices of the controls and external forces. 

3 . Control algorithm 

We seek to minimize the performance index of the form 

T 

J = J(x T QjX + i T Q 2 i + x T Q 3 x + 'i' T Q 4 x + u T Ru)dt (2) 

o 

which is the analogous of the usual linear quadratic cost, where Q,,i = are 
symmetric positive semi-definite matrices and Q 4 and R are symmetric positive definite 
matrices. 

Introducing the change of variable 
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and using the jerk as a control variable via 

, with x = w , system (1') can be transformed into 
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are 3 n x (ti+m) and 3 n x r location matrices respectively for the controls and external 
forces. Here O and I denote respectively the zero matrices and the identity matrix of 
appropriate dimensions. The corresponding performance index is given by 
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This problem can be solved by the Riccati equation approach. We assume that the control 
is regulated by the generalized state vector, i.e., we seek a solution of the form 

a(t) = P(t)z(t) . (13) 

Here P(t) is a symmetric and differentiable matrix. When the external excitation vector is 
neglected, P(t) satisfies the so called Riccati differential equation 
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( 14 ) 


P(t) + P(t)A - — P(t)BR -1 B T P(t) + A T P(t) + 2Q = O, P(T) = O 
2 

The solution is obtained through backward integration in time. Since in most structural 
engineering applications the matrix P(t) remains constant throughout and drops to zero 
near T, P can be assumed constant. It follows that P satisfies the algebraic Riccati 
equation 

PA-^PBR i B t P + A t P + 2Q=o (15) 

u(t) = R~ l B T Pz(t) (16) 

2 

If the solution of the Riccati equation obtained is written the form 
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(17) 


where each P y is an n x n matrix, then 

u(t) = -^R ‘D T M ,X P 2 z(t)= ^R J D T M lT (P 2lXl (t) + P 22 x 2 (t) + P 23 x 3 (t))(18) 
and 

w(t) = Q 4 1 (P 31 x, (t) + P 32 x 2 (t) + P 3 3* 3 (t)) (19) 

hence from (3), the expressions (18) and (19) can be rewritten as 

«(t) = - i R -D* M ‘ x (P 21 x(t) + P 22 x(t) + P^t)) (20) 

and likewise 


w(t) = x(t)= - Iq 4 ‘ (P 31 x(t) + P 32 x(t) + P 33 x(t)) . (21) 

Once the matrix P is obtained and the case the generalized state vector is available for 
measurement the control can be designed by appropriately choosing the weighting 
matrices Q and R 
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4. Applications 

The analysis proposed will lead to the design of more efficient controls for jerky excitation 
like earthquake loads and the design of better controllers to prevent motion sickness. 

5. Conclusion 

The concept of jerk minimization was introduced in this paper with a detail of how this 
promising new concept can be applied to linear systems via quadratic performance index 
involving the cost due to the jerk. This was done through a special change of variable and 
using the jerk as a control variable. The values of the optimal control are obtained using 
the so called of Riccati equation. This result is a fruitful step toward designing efficient 
controllers for the minimization of vibration in structures. 
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